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1. Scientific motivation for studying tropospheric OH 


3. Preliminary results for expected changes in OH TROP 


Sensitivity of OH TROP to climate 

change 

• IPCC (2001) predicted future OH would fall due to rising CH 4 

• Best knowledge of global tropospheric OH (OH TROP ) comes 
from analysis of CH 3 CCI 3 observations 

• The CH 4 lifetime inferred from OH TROP based on CH 3 CCI 3 is 8.9 
years; yet many sources use a CH 4 lifetime of ~12 years, which 
was found by modeling studies described in IPCC (2001) 

• The 12 year CH 4 lifetime is central to subsequent IPCC reports 
and is called the "perturbation lifetime" 

• Our work is motivated by understanding: 

- what factors other than rising CH 4 will affect OH TROP 

- veracity of notion that future OH TROP will decline 


Importance of OH TROP 

• OH serves as the main sink for many species: CO, CH 4 , 

CH 3 CCI 3 , all HCFCs and HFCs, biogenic VOCs (isoprene and 
terpene), anthropogenic VOCs (formaldehyde, benzene, 
toluene, etc.) 

• HO x is integral to 0 3 formation and destruction pathways; 
alone, HO x depletes 0 3 , but with NO x and VOCs, HO x creates 0 3 


Methane Oxidation and OH 


Previous Work on OH TROP Trends 


• Under low NO x , the oxidation of CH 4 destroys HO x 

• Under high NO x , the oxidation of CH 4 produces HO x 

• Loss of methylperoxy radical CH 3 0 2 with NO or H0 2 is critical 

• Future levels of CH 4 and NO x are highly uncertain 

- CH 4 between now and 2100 varies dramatically among 
RCP scenarios 

- NO x between now and 2100 will depend on whether the 
developing world implements selective catalytic reduction 
on coal power plants and catalytic converters on cars 

• Future OH TROP will also depend on overhead column 0 3 , 
local humidity, and biogenic emission of VOCs, etc. 

• Given this complexity, we expect CCMs to project 

a wide range of values for OH TROP during the rest 
of this century 2 ooo 
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• Montzka et o\. ; 2011 found that OH TROP does not vary interannually (from 1997 to present) 

• Prior studies by Prinn et a!., 2001 and Bousquet et a!., 2005 suggest large interannual variability in 
OH TROP (1985-2000) 

• We suggest the OH TROP behavior in all three studies may be physically possible, based on our 
preliminary results 


Current Conclusions from OH TROP Analysis 

• Effect of CH 4 on OH TROP is taken from IPCC (2001), Section 4.2. 1.1, which states that "the feed- 
back of CH 4 on tropospheric OH" found using contemporary chemical transport models is -0.32% 
for every 1% increase in CH 4 (red line, middle panel of figure to right) 

• Primary effect of overhead 0 3 is rise in OH TROP following the 1991 eruption of Mount 
Pinatubo (green dashed line) due to enhanced removal of stratospheric 0 3 by volcanic aerosol 

• Rising H 2 0 from MERRA and AIRS increases OH TROP (orange lines) by an amount comparable 
to the decrease expected from rising CH 4 

• Overall expected change in OH TROP (bottom panel) shows higher level of interannual variability 
prior to ~1999 and lower variability thereafter 
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Globally averaged CH 4 and OH trends vs. Time 
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2. Methods for estimating changes in OH TROP 


4. Future plans for refining estimates of AOH TROP 



Base Case OH Conditions 

• Initial modern-day conditions for OH are taken 
from a time-slice run of GEOS CCM using 2005 
emissions 

• Monthly mean mixing ratios of OH and related 
species are provided on a 144 longitude, 91 
latitude, 72 pressure level grid 

• Calculated changes in OH due to H 2 0 and 
overhead 0 3 are applied to initial OH field 

Annual Average Tropospheric OH column, GEOSCCM 
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MERRA & AIRS H z O 

• Specific humidity files from the NASA Modern-Era 
Retrospective analysis for Research and Applications 
(MERRA) starting prior to 1985, and 


• H 2 0 mixing ratio files from the NASA Atmospheric 
Infrared Sounder (AIRS) starting in 2002 were used to 
establish trends in H 2 0 by latitude 

• OH TROP is assumed to follow the square root of the 
change in H 2 0 using a steady-state assumption 

• We plan to refine the effect of changing H 2 0 on 
OH TROP by examining reaction rates from archived 
runs of GEOS CCM 


Overhead 0 3 Observations 


• Total column 0 3 trends were obtained from the 
NASA merged 0 3 data set, consisting of 
measurements from SBUV, TOMS, and Aura OMI 
instruments 

• We then use our photolysis code to estimate the 
impact on J(0 3 ) -> 0( 1 D) of decreasing initial GEOS 
CCM overhead 0 3 columns by amount suggested by 
the NASA product 

• OH TROP is assumed to change by the square root of 
J(0 3 ) OPD) 



Improve Estimate of 
d(OH TROP ) / d(H 2 0) 

• Reaction rates from recent runs of GEOS CCM 
are archived for reactions such as: 

H 2 0 + O^D) -> 20H 

• Using these reaction rates we will determine 
the proportion of OH that is produced via 
reaction with H 2 0 

• The determined scaling factor would be used to 
calculate a new AOH TROP based on the H 2 0 
trends 

• Estimate time- and pressure-varying values of 
d(OH TROP ) / d(H 2 0) 

• Evaluate discrepancies between MERRA and 
AIRS H 2 0 trends 


Evaluate CH 4 /OH Feedback 

• We will use a box model (details below) to probe 
relationship between CH 4 and OH TROP and its 
dependence on NO x 


Box Model 

We will use the GSFC Combined Stratosphere- 
Troposphere (COMBO) box model provided by 
Chang Lang (JHU): 

• GMI chemical mechanism 

• 118 species, 321 thermal reactions, and 81 
photolysis reactions 

• 5 modules: 

-Aerosol optical depth & surface area 
-Photolysis scheme 
-Thermal reactions scheme 
-Differential eqn solver 
-Input-output 

• Fast-JX photolysis & SMVGEAR II solver 


Propagate Uncertainties 

• Calculate uncertainties in AIRS and MERRA 
H 2 0 and NASA 0 3 product 


• Evaluate standard deviation in average fraction 
of OH production occurring via H 2 0 + O^D) 

• Estimate uncertainty in the box model 
evaluation of d(OH TROP ) / d(CH 4 ) 


Assessing OH TROP in CCMs 

• Through our involvement with the IGAC / SPARC 
Chemistry-Climate Model Initiative, we have 
requested: 

1. hourly, instantaneous output from participating 
CCMs 1 day/season, 1 year/decade 

2. archival of all species, reaction rates, J-values, 
and physical parameters relevant to OH 
chemistry 

3. this output for both the REF-C1 (hindcast) and 
REF-C2 (future) runs 

• We plan to assess the causes of differences 
between OH in the CCMs 

• Use of the box model enables us to distinguish 
between OH differences due to chemical 
mechanism and those due to differences in OH 

precursors 

• We can also use this output to predict future trends 
in OH TROP , based on CH 4 , H 2 0, and overhead 0 3 from 
the future CCM runs 










